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1. Introduction 


QCD could in principle violate CP, as it supports the so-called “theta term” in its Lagrangian: 

^ QCD 9 e 32 ^ ^ Vpa FgvF pa (U) 

which violates CP thanks to symmetry properties of the Levi-Civita epsilon tensor. This term is 
topological, in the sense that its spacetime integral is always an integer multiple of 0, 

2 = 32 tt 2 / d * X ^ Vp %vV G Z (1.2) 

and we call Q the topological charge. The parameter 6 is truly an angle—it can in principle take 
any value only in because the weight in the path integral e lSQCD °< e'@ e is periodic in 0, 

as Q is an integer. From limits on the neutron electric dipole moment, an observable that should 
reflect any strong CP violation, one concludes that |0| < 10 -1 °. The Strong CP Problem asks: if 6 
could naturally take any value in (— n, n], why is it so tiny? 

Axions provide a dynamical mechanism to alleviate this fine tuning. Peccei and Quinn ob¬ 
served that if the parameter 6 were promoted to a dynamical field a with the Lagrangian 
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the Strong CP problem could be solved[3]. The vacuum expectation value of the axion field au¬ 
tomatically adjusts itself so that ( a/f a + 6) vanishes, as can be seen from effective field theory 
arguments. Moreover, the axion’s mass m a satisfies the relation 
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(1.4) 


where Fqcd(T, 0) is the free energy of QCD, which depends on the temperature T and the pa¬ 
rameter 0. The relevant derivative of Fqcd is the temperature-dependent topological susceptibility, 
denoted %(T). Note that m a is not a free parameter of axion physics: given f u and the QCD 
topological susceptibility we know m a , which is temperature dependent. 


2. The Over-Closure Bound 

It was quickly realized that axions not only solve the Strong CP Problem but also are a plau¬ 
sible dark matter candidate[4, 5]. So, it is interesting to use observational knowledge to try and 
constrain axion properties. One such constraint, the over-closure bound, gives a lower bound 
on the axion mass by requiring that axions not account for more than the observed dark matter 
abundance [6]. 

Fligh temperature arguments imply that the QCD topological susceptibility X vanishes at high 
temperatures, so that at the beginning of the universe axions are massless. In particular, at asymp¬ 
totically high temperatures the dilute instanton gas model (DIGM) should be valid, and it predicts 

X{T) - (j/T c y 
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where T c the temperature of the chiral transition can be traded, for example, for Aqcd by simply 
readjusting C, and n is predicted to be exactly 7 when N c = 3. We know from chiral perturbation 
theory (^PT), on the other hand, that at zero temperature the topological susceptibility is[4] 


Z(0) 


m u m d 


(m u + m d ) 


- f 2 m 2 

) J n ,n K 


3.6 • 10 7 MeV 4 


( 2 . 2 ) 


where f K ~ 93MeV, so that % must increase as temperature decreases, but only to a point. 

As cosmological evolution proceeds the universe expands and cools. From the equations of 
motion we can estimate that the axion number becomes fixed when roughly 


3>H(T) k, m a (T) or equivalently 9H 2 {T)f~ « %(T) (2.3) 


which fixes a temperature 7) and a time t\. Put another way, as the universe evolves the topological 
susceptibility becomes large enough that the axion Compton wavelength is comparable to the size 
of the universe, and the axion’s dynamics can kick in. A more complete calculation would solve 
the equations of motion numerically, as is done in Reference [7]. We use 

h 1 { t )= 4 ^- (t)t ' ,2a> 


where Mp\ is the reduced Planck mass \Jhcl%nG ~ 2.4 x 10 18 GeV and g*(T) is the effective 
number of degrees of freedom as a function of temperature given in Ref. [7]. 

After production ceases, axions follow a familiar but slightly modified story — they get di¬ 
luted by the universe’s expansion but, unlike a typical dark matter candidate, continue to get more 
massive until the chiral transition of QCD freezes the topological susceptibility at today’s value 
given in (2.2). An adiabatic analysis shows that[8] 


pW(t? 

m a (t) 


the number of axions in a fixed comoving volume 


(2.5) 


is conserved instead of the usual p(t)R(t) 3 . The cosmological equations of motion provide the 
initial energy density [9], 

p(Tx) = \m 2 a f 2 a el= l -m)6l (2.6) 

which is a function of 6 \—the quantity (a) / f a + 0, the coefficient of the CP-violating term in 
the Lagrangian, at the temperature 7j. If Peccei-Quinn breaking happens after inflation, then we 
can expect different regions of our universe to have different values of 0 1 . Averaging over the 
different regions gives the value (0 2 ) = 7T 2 /3. From now on we will discuss only this late-breaking 
scenario—in a scenario where PQ symmetry breaks early, instead of a lower bound on m a we can 
get a constraint on a combination of m a and the initial condition 6\. 

With mild algebraic massaging, spatially averaging 0], and use of (1.4), we rewrite (2.5) 




(2.7) 


where Ty = 2.73K is the present-day temperature of the CMB and where both factors clearly de¬ 
crease with increasing 7j. We use the ^PT value x(0) as an approximation for x ( Ty ). Cosmology 
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provides R(T), so we now know the present-day axion abundance, given the topological suscepti¬ 
bility x and the temperature T\ (itself given by x, fa, and cosmological inputs). Since (2.3) implies 
that for a fixed topological susceptibility T\ decreases with increasing f a , it is apparent that the 
axion energy density today increases with f a , and that requiring p{Ty) be less than the observed 
dark matter energy density yields an upper bound on f a and subsequently a lower bound, through 
(2.2), on the present-day axion mass. Importantly, this argument also implies that if there is more 
than one kind of dark matter, the lower bound gets tighter. 

Prior state-of-the-art calculations include calculations from the naive dilute instanton gas model 
(DIGM) and numerical calculations from the interacting instanton liquid model (IILM). Refer¬ 
ence [7] found that the IILM predicts m a > 21 ± 2/icV. However, at high temperatures they also 
found that various models differ in their predictions for x{T) by an order of magnitude or more, 
so it is challenging to evaluate the systematic error in their result. To make sure the axion is not 
mistakenly missed, the Axion Dark Matter Experiment (ADMX) takes a very conservative over¬ 
closure bound of ~ 1 jUeV. In principle, a calculation of x directly from lattice QCD would yield 
totally controlled errors and an entirely trustable over-closure bound. 

3. Lattice Calculation 



T/T c T/T c 


Figure 1: The left panel shows a study of discretization effects at three temperatures for all four definitions 
of the topological charge we use. The right panel shows a study of finite volume effects at three different 
temperatures. Each ensemble is slightly offset from its temperature for ease of visibility. The statistical error 
bars on each data point are much smaller than the markers and are suppressed for clarity. 

Unfortunately, a direct calculation from QCD will prove to be extremely expensive. As a pre¬ 
liminary calculation, though, it makes sense to study pure glue—Yang-Mills with no quarks. This 
simplification allows us to take advantage of heatbath algorithm, which enables us to calculate on 
ensembles with huge volumes and gather tremendous statistical samples with short autocorrelation 
times. We studied temperatures in the range 1.2 to 2.5 times T c the critical temperature of pure 
glue (approximately 284 MeV), at two different lattice spacings corresponding to /V-=6 and 8 time 
slices. Previous studies went out to 1.3T C , and our calculation agrees nicely with those results. The 
ensembles are described in detail in Ref. [2], Other groups have since shown similar studies with 
different techniques [10, 11], and while there is some tension in the lattice results, even the most 
conservative (lowest) lower bound excludes all the parameter space that ADMX plans to search 
through 2017. Additional discussion along these lines is postponed until the next section. 
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Figure 2: The left panel shows our best points with a DIGM-inspired (power-law) fit. The right panel shows 
the two sides of (2.3) for three example choices of f a . The intersection point for a given f a sets T\. The 
statistical error bars on each data point are much smaller than the markers and are suppressed for clarity. The 
statistical errors are thinner than the solid line—the band represents a very conservative systematic fitting 
error. More details can be found in Reference [2]. 


The topological susceptibility % can be calculated on a finite volume by measuring the topo¬ 
logical charge Q on each configuration and taking the ensemble average, 


X = lim 

V->°c 


(Q 2 )-(QY 

v 


(3.1) 


We measure Q through the naive discretized analogue of (1.2), 


e. = ^T e " vp Wvn 


pa 


(3.2) 


where □ is the lattice plaquette which reduces to the field strength F in the continuum limit. Be¬ 
cause of the lattice discretization, 2 m is not an integer. We massaged the distribution of Q in four 
different ways (as detailed in Ref. [2]) and the method of Ref. [12] showed very little change as we 
changed the lattice spacing and size, as shown in Figure 1 . Moreover, as the temperature increased 
the other definitions seemed to approach that definition, possibly indicating its reliability. Thus, for 
everything that follows we used the results from that definition. 

We did not study discretization effects at every temperature, so because we observe little lattice 
spacing dependence we simply take the largest volume at the coarser spacing, AT = 6. Those points 
are shown in Figure 2, with a fit inspired by the DIGM in (2.1), where we fit both the constant C and 
the exponent n. At asymptotically high temperatures n = 7, but we find that, at these temperatures 
near T c , n = 5.64 ± 0.04. It is interesting to understand how this difference disappears and %(T) is 
eventually described by a dilute instanton gas[l 1]. 


4. Pure-Glue Bound 

We can use the best DIGM-inspired fit to extrapolate to high temperatures. Of course, eventu¬ 
ally it is important to have the region of cosmological interest under numerical control, especially 
if one desires a rigorous bound from QCD. 
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Plugging %(T) given by our DIGM-inspired fit into (2.3) and (2.7) and requiring that p(Ty) 
accounts for at most the 27% of the universe that is dark matter implies f a < (4.10 ±0.04)10* 1 GcV 
and, through (2.2), the present-day axion mass m a > (I4.6±0.1 j/./cV. This is much higher than the 
bound ADMX uses for guidance, and excludes the masses ADMX plans to search through 2017. 

While our calculation is not a calculation from full QCD, it might still be useful for experi¬ 
mental guidance, as will be explained now. General arguments suggest that x(T) should vanish 
if at least one quark mass vanishes. Pure glue corresponds to the opposite limit, with all quark 
masses going to infinity. It’s extremely plausible, then, that the topological susceptibility for real 
QCD would be below that of pure glue. 1 

We further observe that with either a smaller or faller-fasting topological susceptibility the 
mass bound always gets stronger. We can show this by considering a fixed f a which saturates the 
bound for a given %(T) and asking how p(Ty) changes with a small change in the topological sus- 
ceptiblity x{T) —> x(T) + §x(T). If p(Ty) increases, the previously-allowed /„ will be disallowed 
and the mass bound will be stronger. Equation (2.3) indicates 7) shifts by 57) given by 


52) 

T] 


Sx(Ti) 


4- 


X(L) dT 


& X(Ti) 


(4.1) 


(assuming g*(7)) to be a constant) and thus from (2.7) p(T y ) changes by an amount 5p satisfying 
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(4.2) 


Since the denominator in (4.2) is positive the sign of 5p is controlled by the numerator. In the 
cosmological region of interest the universe is radiation-dominated and thus R o< T 1 , so that 5p 
is positive when 8% is negative—decreasing y^(T) tightens the over-closure bound. 

Since it’s quite plausible that the topological susceptibility of real QCD is below that of pure 
glue, we expect the mass bound from QCD to be a stronger bound than our pure-glue bound of 
m a > (14.6 ± 0.1 )pcV. This strengthens the value of our bound, in that it suggests our bound is 
more conservative than the bound from QCD. 


5. Conclusions and Outlook 

Peccei-Quinn symmetery nicely cleans up the Strong CP problem and simultaneously pro¬ 
vides a plausible, largely unconstrained dark matter candidate—the axion. If the topological sus¬ 
ceptibility of full QCD is indeed lower than that of pure glue, if the estimate (2.3) is reliable, and 
if PQ symmetry broke after inflation, then with its current operational plan ADMX cannot find 
axions until 2018. If other dark matter species account for a sizable portion of the total dark matter 
abundance, the over-closure bound will be tighter and the possibility of experimental observation 
of axions will be further postponed. On the other hand, if an axion is discovered with a mass lower 

'Some speakers[13, 14] at LATTICE 2015 show values of x{T) at 1 < T/Tc < 1.2 for domain wall or overlap 
fermions with light pion masses. 
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than our lower bound ~ 15/ieV, a complete calculation of the over-closure bound from QCD be¬ 
comes urgent, because if that axion should be excluded (as our arguments in the previous section 
suggest) it implies that PQ symmetry broke early—a possible signal of preinflationary dynamics. 

One future direction is the calculation of higher moments of the topological charge at finite 
temperature, which may prove necessary for full numerical control over the cosmological axion 
equations of motion. Additionally, the incorporation of quarks and the measurement of the topo¬ 
logical charge Q via the Wilson flow method are obvious steps towards a fully controlled bound. At 
the temperatures of cosmological interest, topological freezing may make it necessary to measure 
the susceptiblity % using other techniques. Obvious candidates include fixed topology methods and 
open boundary conditions. Finally, QCD with real, finite 6 may be amenable to Langevin methods, 
which would permit a full calculation of the free energy Fqcd(7\ 6). 
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